microcopy  resolution  test  chart 

national  bureau  Of  SIANOAROS  l%i  A 


I 


AD-A191  888 


WHdMITED 


BR 1 04725 


yv-v'-^^iVyr'...:-' 


Report  No.  87013 


IteNALS  AND  RADAR  ESTABLISHMENT, 
MALVERN 


ROYAL  SIGNALS  AND  RADAR  ESTABLISHMENT 
Report  87013 


Title:  The  Avoidance  of  Collisions  for  Newtonian  Bodies 

with  Hidden  Variables 

Author:  B  D  Bramson 

Date:  October  1987 


SUMMARY 


The  collision  avoidance  of  a  pair  of  uniformly  moving  bodies  Is 
considered  In  three  dimensions.  The  relative  motion  of  the  bodies 
yields  an  expression  relating  the  time  to  closest  approach,  the  minimum 
range,  the  current  range  and  Its  rate  of  change,  other  variables  being 
unobservable.  A  Boolean  relation  Is  then  proposed  that  Is  satisfied 
whenever  the  minimum  range  and  time  to  closest  approach  simultaneously 
fall  below  given  thresholds.  The  relation  Is  further  studied.  In 
particular  with  regard  to  the  Issue  of  false  and  premature  alarms. 

An  airborne  collision  avoidance  system  Is  a  possible  application,  f 
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1  INTRODUCTION 


This  paper  concerns  Che  relative  motion  In  three  spatial  dimensions  of 
a  pair  of  uniformly  moving  bodies,  idealised  as  point  particles.  The 
purpose  Is  to  generate  relations  between  observables  that  will  be  useful 
for  the  prediction  and  avoidance  of  collisions.  Applications  Include 
the  design  of  an  airborne  collision  avoidance  system. 

The  problem  would  appear  to  be  mathematically  trivial  save  for  a  crucial 
engineering  constraint.  It  Is  given  that  the  apparatus  for  collision 
avoidance,  carried  by  at  least  one  of  the  bodies.  Is  capable  of 
measuring  directly  only  the  relative  range  of  the  bodies  and  Its  rate 
of  change  In  time.  It  Is  assumed  that  further  time  derivatives  are 
not  computable  with  sufficient  accuracy  and  that  other  observables, 
like  relative  bearing  and  velocity,  are  hidden. 

Section  2  discusses  the  relative  motion  of  the  two  bodies,  expressions 
being  presented  for  the  time,  t,  to  closest  approach,  henceforth  called 
the  warning  time,  and  the  mlnlmtim  range,  m.  Both  m  and  t  depend  on  the 
relative  velocity,  v,  and  are  therefore  unobservable;  but  an  expression 
Is  found,  Independent  of  v,  that  relates  m  and  t  to  the  current  range, 
r,  and  Its  rate  of  change,  r. 

In  section  3.  a  pair  of  threshold  constants  are  Introduced,  mu  and 
tau,  for  m  and  t  respectively,  mu  defines  a  minimum  safe  separation 
for  the  two  bodies  while  tau  defines  a  minimum  warning  time  whenever  mu 
exceeds  m.  Indeed,  a  Boolean  relation,  danger.  Is  defined  between  r, 
f,  m,  t,  mu  and  tau  but,  being  dependent  on  v.  It  Is  unobservable.  An 
expression  Is  further  derived  for  the  time  for  r  to  reach  mu,  given 
appropriate  pre-conditions. 

In  section  4,  an  observable  Boolean  relation,  alarm.  Is  presented 
Involving  r,  r,  mu  and  tau.  Its  strength  lies  In  Its  ability  to  flag 
danger  reliably  while  Its  weakness  lies  In  Its  tendency  to  generate 
either  false  or  premature  alarms. 

Section  5  briefly  Illustrates  the  alarm  In  configuration  space. 


2  UNIFORM  RELATIVE  MOTION 

At  time  u,  consider  a  pair  of  uniformly  moving  point  particles  with 
relative  position  vector  X(u)  and  (constant)  relative  velocity  vector 
X.  (’  denotes  d/du.)  Let  their  relative  range  be  r(u) ,  being  the 
magnitude  of  X(u),  so  that  their  relative  range  rate  Is  r(u).  Further, 
let  their  relative  speed  be  v,  being  the  magnitude  of  X. 

To  calculate  the  warning  time,  t(u),  and  the  minimum  range,  m,  the 
equation  of  relative  motion  Is  used.  This  determines  the  relative 
position  vector  at  time  u  from  that  at  time  0  according  to 

X(u)  -  X(0)  +  X  u  ,  (2.1) 

whence  the  behaviour  of  the  relative  range  Is  given  by 

r*(u)  -  r'(0)  +  2r(0)  i(0)  u  +  v*  u*  .  (2.2) 

Henceforth  assume  that  v  /-  0.  This  allows  us  to  choose  the  origin  of 
time  to  correspond  to  the  Instant  of  closest  approach.  This  being  so. 


r(0)  -  m,  r(0)  -  0  and  t  -  -u.  Thus,  dropping  the  argument  (u) , 

r*  “  +  V*  u*  .  (2.3) 

Differentiating  with  respect  to  u  and  using  t  -  -u  provides  an 
expression  for  the  warning  time, 

t  -  -  r  i  /  V*  .  (2.4) 

while  the  minimum  range,  m,  is  given  by 

m*  -  r*  (1  -  r*  /  V*)  .  (2.5) 

Note  from  (2.4)  that 

(r  >  0  AND  r  <  0)  -  t  >  0  .  (2.6) 

For  V  “  0,  m  and  t  are  both  undefined.  Even  for  v  /-  0,  neither  m 
nor  t  are  known  since  v  is  a  hidden  variable.  Nevertheless,  they 
possess  a  relation  that  is  independent  of  v  and  this  follows  from 
equations  (2.3)  and  (2.4): 

J  2  • 

m  “  r  +  r  r  t  .  (2.7) 

Vhen  the  bodies  are  converging  in  range,  r  <  0,  it  follows  that 

m*  +  r  |r|  t  »  r^  AND  t  >  0  (2.8) 

which,  for  given  r  and  r  (<  0) ,  is  Illustrated  in  figure  1. 


t  " 


■  —  ■  ■  ■  ■  ■  > 

r  m 

Figure  1:  Knowing  r  and  r  (<  0)  yields  a  relation  between 
the  warning  time  t  and  the  minimum  range  m. 

3  THE  IMPOSITION  OF  THRESHOLDS 

Let  us  now  impose  a  minimum  safe  separation  on  the  two  bodies,  call  it 
mu  (>  0) .  We  require  that  r  >  mu  for  all  time.  Further,  if  m  <-  mu  <  r 
and  r  <  0,  we  shall  wish  to  be  alerted  should  the  warning  time  be  less 
than  a  certain  threshold,  tau  (>  0),  this  being  considered  dangerous. 
Within  this  context  it  is  vital  to  have  a  rigorous  definition  of 
danger.  So,  consider  the  time  dependent  Boolean  expression 
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danger  - 

(r  <-  mu  OR  r  <  0  AND  m  <-  mu  AND  t  <-  tau)  . 


(3.1) 


(AND  binds  tighter  than  OR.) 

Sadly,  Tinless  r  <-  mu,  danger  Is  unobservable  being  dependent  on  v 
and  an  attempt  to  resolve  this  Is  made  In  the  following  section. 

To  conclude  this  section,  we  calculate  the  time,  t_mu,  for  r  to  reach 
the  critical  value,  mu,  given  that  m  <-  mu. 

Using  equation  (2.3)  we  see  that  If  m  <-  mu,  r  reaches  mu  when 

u*  -  (mu*  -  m*)  /  V*  .  (3.2) 

The  two  solutions  for  m  <  mu  correspond  to  events  before  and  after 
closest  approach.  It  follows  that  the  time  to  reach  threshold  before 
the  Instant  of  closest  approach  Is  given  by 

t_mu  -  t  -  (sqrt(mu*  -  m*))  /  v  .  (3.3) 

The  relative  motion  of  the  two  bodies  Is  Illustrated  In  figure  2. 


<  . V  t . > 

<  - V  t  mu - > 


Figure  2:  The  motion  of  1  relative  to  2  when 
r  <  0  AND  m  <«  mu  <  r. 


4  AN  OBSERVABLE  ALARM 

The  Boolean  relation,  danger,  defined  In  equation  (3.1)  Is  unobservable; 
but  In  section  2  we  derived  a  relation  (2.7)  between  m  and  t  In  terms 
only  of  observable  quantities.  Constraining  m  and  t  as  In  (3.1)  Imposes 
constraints  on  r*  from  (2.7).  Accordingly,  define  the  Boolean 
expression,  alarm,  by 

alarm  -  ^  ^  ^  (4.1) 

(r  <-  mu  OR  r  <  0  AND  r*  <-  mu  -  r  r  tau)  . 

Using  the  relations  (3.1),  (4.1)  and  (2.8),  It  follows  that 

danger  — >  alarm  .  (4.2) 

Unlike  danger,  alarm  Is  observable  and  the  question  therefore  arises 
as  to  whether  It  may  be  used  In  place  of  danger.  To  examine  this 
possibility,  we  use  equations  (2.3),  (2.4)  and  (2.6)  together  with  the 
fact  that  u  -  -t,  to  rewrite  (4.1)  as 


irHIWKIUI.WWM  ■W  V  Nt  m  ■> 


alarm  -  (^-3) 

(m‘  +  v^  t*  <-  mu  OR 
t  >  0  AND  m*  +  ▼*  t(t  -  tau)  <-  mu  )  . 

To  examine  this  further,  take  the  conjunction  with 

t  <-  0  OR  t  >  0  (A. 4) 


and  use  the  property  that 

(t  >  0  AND  m*  +  V*  t*  <-  mu^  )  — > 

(t  >  0  AND  m*  +  V*  t(t  -  tau)  <-  mu*)  . 


This  yields  the  expression 


alarm  - 

(t  <-  0  AND  m*  +  V 
t  >  0  AND  m*  +  V* 
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t*  <-  mu^  OR 
(t  -  tau/2)*  <-  mu* 


+  V 


tau* /A) 


(4.5) 


(4.6) 


which  defines  the  union  of  the  interiors  of  quadrants  of  two  ellipses 
in  the  (t,  m)  plane,  with  centres  (0,  0)  and  (tau/2,  0)  respectively, 
and  is  Illustrated  in  figure  3. 


Figure  3:  The  shaded  regions  comprise  alarm.  The  arrow 
of  time  is  the  direction  of  increasing  u. 


Note  that  the  alarm  is  raised  only  if  m  is  bounded  above  by  a  certain 
▼elocity  dependent  expression: 

alarm  — >  m  <■•  max  ,  (4.7) 

with  max  ->  sqrt(mu*  +  v*  tau  /4)  .  (4.8) 

Further,  the  boundary  of  alarm  includes  the  points: 

A  (0,  mu)  ,  (4.9) 

B  (tau/2,  max) 

C  (tau,  mu)  , 

D  (tau/2  ♦  max/v,  0) 

F  (-mu/v,  0) 


and 


In  fact  alarm  comprises  Che  disjunction  of  four  Boolean  expressions 
corresponding  to  the  four  distinct  regions  In  figure  3.  These  may  be 
Identified  formally  by  forming  the  conjunction  of  (4.3)  with  the 
expression 

(m  <-  mu  OR  m  >  mu)  AND  (4.10) 

(t  <-  0  OR  0  <  t  <-  tau  OR  t  >  tau)  , 


yielding 


alarm  -  ,  ,  .  ,  (4.11) 

1  (t  <<■  0  AND  m  +  V  t  <-  mu*  OR 

2  0  <  t  <-  tau  AND  m  <>  mu  OR 

3  0  <  t  <  tau  AND  m  >  mu  AND  m^  +  v*  t(t  -  tau)  <-  mu*  OR 

4  t  >  tau  AND  m  <  mu  AND  m*  +  v*  t(t  -  tau)  <-  mu*)  . 

Region  1  Is  bounded  by  the  lines  OA  and  OF  and  Che  ellipse  segment  AF 
while  region  2  Is  bounded  by  the  rectangle  OACE.  The  union  of  regions 
1  and  2  correspond  to  danger.  From  (4.2)  we  expect  the  alarm  to  be 
raised  here.  The  remainder  of  alarm  comprise*;:  region  3,  a  false  alarm, 
bounded  by  the  ellipse  segment  ABC  and  the  line  AC;  and  region  4,  a 
premature  alarm,  bounded  by  the  ellipse  segment  CD  and  by  the  lines  CE 
and  ED. 
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We  now  examine  the  behaviour  of  alarm  for  the  two  cases ,  m  <-  mu  and 
m  >  mu.  In  each  case  consider  the  motion  of  a  point  In  the  (c,  m) 
plane,  parallel  to  the  t  axis  and  In  the  direction  of  t  decreasing, 
representing  the  relative  motion  of  the  two  bodies. 

CASE  1  m  <«  mu.  The  alarm  comes  on  when  the  warning  time  satisfies 

t  -  tau/2  +  1/v  sqrt(max^  -  m*)  >«■  tau  ,  (4.12) 

corresponding  to  a  point  on  the  ellipse  segment  CD,  and  goes  off  when 

t  «■  -1/v  sqrt(mu*  -  m*)  ,  (4.13) 

corresponding  to  a  point  on  the  ellipse  segment  AF. 

When  m  -  mu  the  motion  follows  CA,  the  alarm  coming  on  when  t  -  tau  and 
going  off  when  t  <■  0, 
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When  m  <  mu  a  premature  warning  Is  given  when  the  motion  encounters  the 
segment  CD.  Indeed,  the  alarm  comes  on  for  a  value  of  t,  given  by 
equation  (4.12),  that  may  be  arbitrarily  large  for  v  small  enough, 
approximating  to  mu/v  for  a  slow,  collision  course  (along  DEOF) . 

However,  the  time  t_mu  for  r  to  reach  Its  threshold  mu  Is  bounded  above 
by  tau.  To  see  this,  use  equations  (3.3)  and  (4.3)  to  show  that,  on  the 
segment  CD, 


s.' 

s. 

A. 


t_mu  -  t  -  sqrt(t(t  -  tau))  AND  t  >-  tau 


(4.13) 


t_mu  decreases  monotonlcally  from  tau  as  t  Increases  from  tau  and  Is 
bounded  below  by  tau/2  (see  appendix).  The  expression  (4.13)  Includes 
the  case  when  m  «  mu  (t_mu  -  t  ->  tau)  and  Is  Illustrated  In  figure  4. 


After  region  4  the  motion  proceeds  Into  region  2  and  thence,  following 
the  Instant  of  closest  approach.  Into  region  1  where  the  bodies  while 
diverging  In  range  are  still  within  a  distance  mu  apart. 


Figure  4:  When  m  <•>  mu.  the  elarm  is  raised  as  soon  as  t_mu 
reaches  a  value  that  satisfies  tau/2  <  t_mu  <-  tau. 


CASE  2  ffl  >  mu.  This  appears  to  be  a  false  alarm.  The  alarm  Is  raised 
If  m  <-  max,  defined  In  (4.8).  It  comes  on  when 

t  ■  tau/2  +  1/v  8qrt(max^  -  m*)  ,  (4.14) 

which  satisfies  tau/2  <>  t  <  tau,  and  goes  off  when 

t  -  tau/2  -  1/v  8qrt(max^  -  m*)  ,  (4.15) 

satisfying  0  <  t  <  tau/2. 

For  V  small  compared  with  mu/tau,  the  upper  bound  for  m,  max,  achieved 
at  B,  lies  close  to  mu  while  for  large  v  it  approximates  to  v  tau/2. 

To  conclude  this  section  a  brief  comparison  will  be  made  between  alarm 
and  one  version  of  the  modified  tau  criterion  [1].  Versions  differ  on 
their  behaviour  following  the  Instant  of  closest  approach,  the  one 
chosen  here  being  given  by 

modlfled^tau  -  (4.16) 

(r  <-  mu  OR  r  <  0  AND  r  <-  mu  -  r  tau)  . 


THEOREM  alarm  — >  modlfled_tau  . 

In  other  words,  modlfled_tau  Is  raised  whenever  alarm  Is  raised, 
follows  swiftly  from  the  definition  (4.1)  on  noting  that 

z 

r  >  mu  — ■>  mu  -  r  r  tau  <  r  (mu  -  r  tau)  . 


(4.17) 

This 


(4.18) 


5  CONFIGURATION  SPACE  PICTURE 

The  analysis  in  the  previous  section  concerned  relations  between  the 
minimum  range,  m,  the  warning  time,  t,  and  the  time,  t_mu,  for  the  range 
threshold  to  be  reached,  when  the  alarm  is  raised.  Following  the 
spirit  of  work  by  Ford  [1,  2],  it  is  also  convenient  to  Illustrate  the 
boundary  of  alarm  In  configuration  space. 


Using  X  and  X  to  denote  the  relative  position  and  velocity  vectors  at 
time  u,  as  in  section  2,  and  defining  max  by  (4.8),  the  expression 
(4.1)  for  alarm  may  be  rewritten 

alarm  -  (|XI  <«  mu  OR  .  (5.1) 

X.x  <  0  AND  |X  -t-  tau/2  X|  <-  max)  ; 

which  represents  the  union  of  the  interiors  of  a  pair  of  spheres  whose 
radii  are  mu  and  sqrt(fflu*  v^  tau*/4)  and  whose  centres  lie  in  the 
direction  of  relative  motion,  a  distance  v  tau/2  apart.  This  is 
illustrated  in  figure  5  from  which  many  of  the  properties  previously 
derived  receive  a  ready  interpretation. 


of  a  pair  of  spheres  in  configuration  space. 


6  CONCLUDING  REMARKS 

This  paper  has  discussed  the  collision  avoidance  of  two  Newtonian  bodies 
in  uniform  motion  when  only  the  relative  range  and  its  rate  of  change 
are  observable.  In  section  3  the  concept  of  danger  was  defined  in  terms 
of  given  threshold  constants,  mu  and  tau,  by 

danger  -  (range  <-  mu  OR  (6.1) 

range  rate  <  0  AND  minimum  range  <-  mu 

AND  warning  time  <-  tau)  . 

However,  danger  is  in  general  unobservable,  the  minimum  range,  m,  and 
warning  time,  t,  depending  on  the  relative  velocity,  v.  Hence,  an 
observable  Boolean  relation,  alarm,  was  Introduced  in  section  4: 

alarm  -  (r  <-  mu  OR  (6.2) 

r  <  0  AND  r*  <“  mu^  -  r  r  tau)  ; 

and  this  is  raised  whenever  there  is  danger: 

danger  •— >  alarm  .  (6.3) 

The  question  then  arose  as  to  the  use  of  alarm  to  indicate  danger. 

In  particular,  if  the  two  bodies  are  converging  in  range  (r  <  0)  what 
are  the  implications  of  alarm?  The  conclusions.  Illustrated  in  figures 
3  and  4,  amount  to  two  Boolean  cases  in  addition  to  danger.  One  of 
these  involves  the  time  t  mu  for  r  to  reach  its  threshold  mu. 
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To  summarise,  for  m  <••  mu  the  alarm  comes  on  with  t  >-  tau  as  soon  as 
t_Bu  becomes  small  enough.  Its  value  lying  between  tau/2  and  tau  and 
depending  on  v;  and  this  will  happen  before  danger  Is  reached.  Whether 
this  Is  regarded  as  a  premature  alarm  depends  on  whether  a  warning  that 
r  Is  approaching  mu  Is  considered  useful.  After  the  Instant  of  closest 
approach  the  alarm  goes  off  as  soon  as  r  exceeds  mu. 

X  Z  2. 

Secondly,  an  alarm  Is  raised  for  mu  <  m  <••  sqrt(mu  +  v  tau  /4) 
as  soon  as  t  becomes  small  enough.  Its  value,  dependent  on  v,  lying 
between  tau/2  end  tau.  The  expression  bounding  m  here  behaves  like  mu 
for  V  small  compared  with  mu/tau  and  like  v  tau/2  for  large  v.  Finally, 
the  alarm  goes  off  for  some  value  of  t  lying  between  0  and  tau/2 . 
Whether  this  case  Is  regarded  as  a  false  alarm  depends  on  whether  a 
level  of  protection  that  Increases  with  v  Is  considered  useful. 

With  regard  to  Its  use  for  a  collision  avoidance  system.  Ford  [2]  has 
compared  the  behaviour  of  alarm  with  that  of  other  criteria.  In 
relation  to  this  It  was  proved,  at  the  end  of  section  4,  that 

alarm  — >  modlfled_tau  ,  (6.4) 

where  modifled_tau  -  (6.5) 

(r  <“  mu  OR  r  <  0  AND  r  <-  mu  -  r  tau)  . 

Thus  modlfled_tau  Is  raised  whenever  alarm  Is  raised.  However,  the 
converse  Is  not  true;  modlfled_tau  Is  raised  more  frequently  than  alarm. 
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APPENDIX 

Consider  the  behaviour  of  the  real-valued  partial  function 

f(t)  «  IF  t  >-  tau  THEN  t  -  sqrt(t(t  -  tau))  ,  (Al) 

with  tau  (>  0)  regarded  as  fixed.  To  examine  f,  put 

f (t)  "  g(*)  with  t  -  tau  cosh*  x  ,  for  x  >-  0  .  (A2) 

As  X  Increases  from  0,  so  t  Increases  from  tau.  Standard  properties  of 
hyperbolic  functions  may  then  be  used  to  show  that 

g  -  tau/2  (1  +  exp(-2x))  ,  (A3) 

which  decreases  monotonlcally  from  tau,  tending  to  tau/2. 


OOCUUrNT  CONTROL  SHEET 


OveriTT  security  class! Ti cation  ot  sheet 


UNCLASSIFIED 


(As  far  as  possible  this  sheet  should  contain  onTy  unclassified  inforaaiion.  If  it  is  necessary  to  enier 
classified  inforaation,  the  bo>  concerned  oust  be  parked  to  indicate  the  classification  eg  (R)  (C)  or  (S)  ) 


1.  DfilC  Reference  (if  knoun)  2.  Originator's  Reference  3.  Agency  Reference  Report  Security 

Report  87013  U/C  Classification 


5.  Originator's  Code  (if  6.  Originator  (Corporate  Author)  Naae  and  Location 

known) 

7784000  RSRE  St  Andrews  Road,  Malvern,  Worcs.  WR14  3PS 


5a.  Sponsoring  Agency's 
Code  (if  known) 


6a.  Sponsoring  Agency  (Contract  Authority)  Naae  and  Location 


7.  Title 


THE  AVOIDANCE  OF  COLLISIONS  FOR  NEWTONIAN  BODIES  WITH  HIDDEN  VARIABLES 


7a.  Title  in  Foreign  Language  (in  the  case  of  translations) 


7b.  Presented  at  (for  conference  napers)  Title,  place  and  date  of  conference 


6.  Author  1  Surnaae,  initials!  9(a)  Author  2 


Bramson,  B.D. 


11.  Contract  Nuaber 


12,  Period 


9(b)  Authors  3. A. 


13.  Project 


10.  Date  pp.  ref. 

10.1987  8 


1A.  Other  Reference 


15.  Distribution  stateaent 


Unlimited  Distribution 


Descriptors  (or  keywords) 


continue  on  separate  piece  of  paper 


Ab»tract  collision  avoidance  of  a  pair  of  uniformly  moving  bodies  is 

considered  in  three  dimensions.  The  relative  motion  of  the  bodies  yields 
an  expression  relating  the  time  to  closest  approach,  the  minimum  range,  the 
current  range  and  its  rate  of  change,  other  variables  being  unobservable. 

A  Boolean  relation  is  then  proposed  that  is  satisfied  whenever  the  minimum 
range  and  time  to  closest  approach  simultaneously  fall  below  given  thresholds. 
The  relation  is  further  studied,  in  particular  with  regard  to  the  issue  of 
false  and  premature  alarms.  An  airborne  collision  avoidance  system  is  a 
possible  application. 
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